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ABSTRACT 


Digital  computer  programs  are  given  for  the  transmission  and 
reflection  coefficients  of  low-loss,  multilayer,  radome  sandwiches. 
Using  matrix  multiplication,  these  programs  yield  design  data  for 
parallel  and  perpendicular  polarization  in  the  most  efficient  manner. 
Typical  data  are  included  for  a  five-layer  sandwich  for  normal  and 
oblique  incidence  at  frequencies  from  zero  to  40  gigacycles. 

Matrix  multiplication  is  also  employed  to  calculate  the  phase 
velocity  of  the  surface  waves  which  may  propagate  on  a  plane  multilayer 
sandwich.  Typical  results  are  given  for  three-layer  and  five-layer 
sandwiches. 

l 

l 

The  transmission  coefficient  of  an  inhomogeneous  plane  layer  is 
determined  with  step-by-step  numerical  integration  of  the  wave  equation. 
A  computer  program  is  included  for  this  problem,  along  with  typical 
computed  data  for  an  exponentially  inhomogeneous  layer. 
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CALCULATION  OF  TRANSMISSION  AND  SURFACE-WAVE 
DATA  FOR  PLANE  MULTILAYERS  AND 
INHOMOGENEOUS  PLANE  LAYERS 

I.  INTRODUCTION 

To  facilitate  the  investigation  of  electromagnetic  windows  and 
radomes,  digital  computer  programs  have  been  developed  for  the  trans¬ 
mission  and  reflection  coefficients  of  low-loss,  multilayer,  radome 
sandwiches.  Using  matrix  multiplication,  these  programs  yield  design 
data  for  parallel  and  perpendicular  polarization  in  the  most  rapid  and 
efficient  manner.  They  have  been  employed  to  study  the  performance  of 
a  five-layer  unsymmetrical  sandwich  for  normal  and  oblique  incidence 
over  a  broad  band  of  frequencies. 

An  understanding  of  the  electromagnetic  effects  of  dielectric 
windows  requires  some  knowledge  of  the  surface  waves  that  may  be 
excited  in  the  windows  by  a  radar  antenna.  This  knowledge  can  be  used 
to  minimize  the  excitation  of  such  waves  and  to  predict  the  effects  of  the 
modes  that  are  excited.  Therefore,  a  computer  program  was  developed 
for  identifying  all  the  modes  that  can  propagate  on  a  given  multilayer 
sandwich,  and  for  calculating  the  phase  velocity  of  each  propagating  mode. 
The  application  of  the  matrix  formulation  to  surface -wave  analysis  is 
believed  to  be  new. 

Another  computer  program  produces  transmission-coefficient  data 
for  inhomogeneous  radomes  in  which  the  dielectric  constant  varies  con¬ 
tinuously  as  a  function  of  distance  from  the  radome  surface.  Discon¬ 
tinuities  are  allowed  for  at  both  surfaces.  Since  the  program  is  based 
on  step-by-step  numerical  integration,  accurate  data  can  be  obtained 
even  for  rapidly  varying  dielectric  constants. 

II.  TRANSMISSION  COEFFICIENTS  OF  LOW- LOSS 
PLANE  MULTILAYERS 

Although  the  matrix  formulation  for  the  analysis  of  plane  multilayers 
is  well  known,  ’  ’  an  outline  of  the  theory  is  included  here  to  facilitate 
the  definition  of  the  quantities  in  the  computer  programs. 
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Consider  a  harmonic  plane  wave  obliquely  incident  on  the  surface  of 
a  stack  of  plane,  homogeneous,  isotropic,  dielectric  slabs  of  finite 
thickness  and  infinite  width  as  shown  in  Fig.  1.  Assume  the  surrounding 
medium  is  free  space. 


Fig.  1.  Plane  wave  incident  on  plane  multilayer . 


In  the  case  illustrated,  the  wave  is  said  to  have  "perpendicular 
polarization"  since  the  electric  field  intensity  vector  is  perpendicular 
to  the  "plane  of  incidence"  (defined  by  the  propagation  axis  of  the  in¬ 
cident  wave  and  the  normal  to  the  dielectric  surface).  Let  the  symbols 
Ej  and  Er  represent  the  electric  field  intensities  of  the  incident  and 
reflected  plane  waves  evaluated  at  the  "incident  point,  "  point  P  in 
Fig.  1;  and  let  E^.  represent  the  electric  field  intensity  of  the  transmitted 
plane  wave  at  the  "normal  exit  point"  Q.  The  reflection  coefficient  R 
and  the  "normal  transmission  coefficient"  Tn  of  the  plane  multilayer 
are  defined  by 


(perpendicular  polarization); 
(perpendicular  polarization). 


(1) 


(2) 


R 


Tn  = 


Er(P) 

Ei(P) 

Et(Q) 

Ej(P) 
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It  is  also  useful  to  define  the  "insertion  transmission  coefficient" 
T  as  follows: 


(3) 


(perpendicular  polarization) 


jkd  cos  6 

1  n  e 


where  d  is  the  total  thickness,  0  is  the  angle  of  incidence  measured 
from  the  normal,  and  k  is  the  free- space  phase  constant 

The  resultant  field  in  each  layer  consists  of  an  outgoing  plane  wave 
and  a  reflected  plane  wave.  Let  the  complex  constants  An  and  Cn  represent 
the  electric  field  intensity  Ex  of  the  outgoing  wave  in  layer  n,  evaluated 
at  its  two  boundaries.  Similarly,  let  Bn  and  Dn  represent  the  reflected 
electric  field  intensity  at  the  two  boundaries  in  layer  n  as  shown  in  Fig.  2. 
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Space 
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Fig.  2.  Plane  wave  incident  on  plane  multilayer . 
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The  field  in  an  arbitrary  layer,  say  layer  n,  is  given  by 

(4)  Ex  =  (ae“YnX  +  be7nx)  e-JkY  sin  0  ^ 

The  propagation  constant  Yn  is  expressed  in  terms  of  the  attenuation 
constant  On  and  the  phase  constant  as  follows: 

I5'  >„  =  “»  +  JPn  • 

Let  us  assume  the  permeability  of  each  layer  is  real  and  the  complex 
permittivity  is  expressed  as  follows: 


(6)  e  =  e'(l  -  j  tan  5). 

From  the  wave  equations  and  Eqs.  (4),  (5),  and  (6)  it  is  found  that 


a  =  (k/J- 2~) |  (prer  -  sin29)2  +  ([i^tan  6)2  -  (jire^  -  sin20) 


and 

(8)  P  =  (k / J~2)  J(jire,r  -  sin2 Q)2  +  tan  &)2  +  -  sin20); 

where  |ir  and  £r  are  the  relative  permeability  and  permittivity: 


(9)  txr-|i/p0 


and 


By  evaluating  Ex  in  Eq.  (4)  at  the  two  points  (x,y,z)  =  (0,0,Zj)  and 
(0,0,  z2),  where  Zj  and  z2  are  the  coordinates  of  the  left  and  right  boundaries 
of  layer  n,  it  can  be  shown  that 

(H)  An=Cne'Yndn 
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and 

(12)  Bn  =  Dn  eVn 


where  dn  is  the  thickness  of  layer  n.  Equations  (11)  and  (12)  can  be 
expressed  by  the  following  matrix  equation: 


Let  t.i  „  and  r  ,,  denote  the  interface  transmission  and  re- 
n+l,n  n+l,n 

flection  coefficients  for  a  wave  in  layer  n+1  incident  on  the  boundary  of 
layer  n.  Further,  let  tn  n.|q  and  rn  represent  the  interface  coeffi¬ 
cients  for  a  wave  in  layer  n  incident  on  the  boundary  of  layer  n+1.  In 
terms  of  these  coefficients,  the  electric  field  intensities  of  the  outgoing 
and  reflected  waves,  evaluated  at  both  sides  of  this  boundary  (between 
layers  n  and  n+1),  are  related  linearly  as  follows: 


(14) 


Cn  =  tn+l,nAn+l+rn,n+l  Dn 


and 


(15) 


Bn+1  ^.n+l^n  +  rn+l,nAn+l  ' 


These  relations  follow  from  the  superposition  theorem  and  the  definitions 
of  the  interface  coefficients. 

It  can  be  shown  that 


<16>  rn,n+l  =  "  rn+l,n  ’ 

<17>  tn+l,n  =  1  +  rn+l.n  ’ 
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(18) 


*"n,  n+1  ^  rn  >  n+1  ^  r  a+1  >  n 

and 

In+l.n  tn, n+1  "  rn+l,n  rn,n;l  ”  1 


By  using  Eqs.  (16)  through  (19)>  Eqs.  (14)  and  (15)  c~n  be  rearranged 
as  follows: 


(20) 

and 

(21) 


^n  “  (A n+1  rn,n+l  ^n+l^n.n+l 


Dn  ~  (Bn+1  '  rn+l,nAn+l)/tn,n+l 


These  can  be  expressed  in  matrix  form  as  follows: 


(22) 


Jn  \  1 


DnJ  ^Tl  \-rn+1>n  1 


"  rn+l,  n /An+1 


B 


n+L 


The  matrix  Eqs.  (13)  and  (22)  can  be  combined  to  obtain  the 
following: 


(23) 


Cn-A  1 


;"Yndn 


'  rn, n-1 e 


y  d 


n  n 


Dn_/tn-l,nUn,n-le"Yndn  eYndn 


D 


Let  us  denote  the  two-by-two  matrix  in  Eq.  (23)  by  the  symbol  MR: 


(24) 


-v  dn 
/  0  ln  n 


•  r. 


a,  n-1 


y  d  . 
„  *n  n 


ie‘Yndn  Yndn 
n,  n-1  e  e 


Repeated  application  of  Eq.  (23)  yields  the  following  matrix  relation 
among  the  electric  field  intensities  A^.^  and  B^j  at  the  "incidence 
surface"  z  =  0,  and  the  intensities  C0  and  D0  at  the  "exit  surface"  z=d:’‘c 

!,!The  dots  in  Eq.  (25)  denote  matrix  multiplication. 
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(25) 


fc  o\ 

(  AN+l^ 

(  =  (1/t)  M.  •  M,  •  M3  •  M.,  •  • 

•  mn  • s  •  | 

\D  0/ 

y  bn+i  j 

where  N  represents  the  total  number  of  layers,  S  denotes  the  matrix 

/i 


(26) 


S  = 


-  r 


N+l.N 


'  ~rN+l,  N  1 


and 

t  =  *0,1  ti,  2  42,3  '  “  tN,  N+l 


Equation  (25)  applies  even  in  the  general  case  where  plane  waves 
are  incident  on  both  outer  surfaces  of  the  multilayer,  provided  the  two 
incident  plane  waves  have  the  same  frequency,  angle  of  incidence,  and 
polarization.  In  the  situation  used  to  define  the  transmission  and  re¬ 
flection  coefficients  of  the  structure,  a  wave  of  unit  amplitude  is  assumed 
to  be  incident  on  one  outer  surface  only  (say  surface  z  =  0),  so  that 


(28)  ^N+l 

(39)  Bn+1  =  R, 

(30)  C0  =  T„, 

and 

(31)  D0  =  0. 

Thus  Eq.  (25)  becomes 


T  \ 

n]  =  (1/t)  • 

0  ' 


m2 


M3  *  M4  •  •  *■  Mj\j  *  S 


(32) 


The  Solution  for  "parallel  polarization"  (where  the  electric  field 
intensity  vector  is  parallel  with  the  plane  of  incidence)  is  obtained  by 
applying  the  theorem  of  duality  to  the  solution  given  above.  Thus,  the 
reflection  and  transmission  coefficients  are  defined  by 


(33) 


-  Hr(P) 


R  =  iix 


Hj(P) 


(parallel  polarization) 


and 

(34) 


_  Ht(Q) 
"  Hj(P) 


(parallel  polarization). 


The  matrix  equations  given  above  apply  also  for  parallel  polarization, 
in  which  case  the  complex  constants  An,  Bn,  Cn,  and  Dn  represent  the 
amplitudes  6  the  magnetic  field  intensities  Hj.  of  the  traveling  waves  in 
layer  n.  Equations  (16)  through  (19)  also  apply  for  parallel  polarization, 
in  which  case  the  interface  reflection  and  transmission  coefficients  are 
defined  by  tfje  ratios  of  the  magnetic  field  intensities  H^.  The  interface 
reflection  coefficients  are  given  by3 


(35) 


]  n+1,  n 


_  ^n  ^n+1  "  4n+l 
^n  ^n+l  *  ^n+l ^n 


(perpendicular  polarization) 


and 


(36) 


r  . .  =  en^n+l  "  £n+l^n 


in+l ,  n 


e'n^n+l  +  ^n+1  Yn 


(parallel  polarization) 


where  V  is  given  by  Eqs.  (5),  (7),  and  (8)  if  the  permeability  of  each 
layer  is  real. 


After  the  indicated  matrix  multiplications  in  Eq.  (32)  are  performed, 
and  the  division  by  t  has  been  carried  out,  the  equation  has  the  form 


(37) 
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Thus 


(38)  Tn  =  a  +  bR 
and 

(39)  0  =  c  +  dR  . 


In  view  of  the  reciprocity  theorem,  the  transmission  coefficient  of  a 
given  multilayer  is  unaffected  by  a  reversal  in  the  direction  of  propagation 
of  the  incident  plane  wave.  If  the  incident  wave  direction  is  reversed 
(see  Fig.  1  or  Fig.  2),  Eqs.  (28)  through  (31)  are  replaced  by  the  following: 

(40)  C0  =  R7  , 


(41) 

(42) 
and 

(43) 


D0  -  i> 
an+i  “  °* 

BN+1  "  Tn 


where  R  denotes  the  reflection  coefficient  for  a  wave  in  free  space 
incident  on  the  right-hand  surface  of  the  multilayer  in  Fig.  2.  (In 
general  R/  differs  from  R,  where  R  is  the  reflection  coefficient  for  a 
wave  in  free  space  incident  on  the  left-hand  surface  of  the  multilayer.  ) 
From  Eqs.  (25)  and  (40)  through  (43), 


where  a,  b,  c,  and  d  are  the  same  as  in  Eq.  (37).  Thus, 

(45)  R  =  bTn, 

(46)  1  =  dTn  , 
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(47) 

ad  -  be  =  1, 

(48) 

Tn  =  1/d, 

(49) 

R  =  -c/d, 

and 

(50) 

R  =  b/d  . 

In  comparison  with  other  methods,  the  repetitive  matrix  multi¬ 
plication  indicated  in  Eq.  (32)  is  systematic  and  easy  to  program  on  a 
digital  computer.  It  is  efficient  with  respect  to  memory  storage  re¬ 
quirements  and  running  time.  Figures  3  and  4  show  Fortran  programs 
(OSU  Version  2)  for  the  transmission  and  reflection  coefficients  of  low- 
loss  plane  multilayers  for  perpendicular  and  parallel  polarization,  using 
the  matrix  formulation.  These  computer  programs  apply  to  non-magnetic 
media  for  which  pL  =  jiQ  and  pr  =  1.  To  reduce  the  memory  storage  re¬ 
quirements  and  the  running  time,  the  programs  utilize  the  following 
approximations  for  low-loss  media  such  as  those  used  in  radome  con¬ 
struction. 


(51) 


rn+l,  n 


J£n+1  -  sinZ9  -  J£n  -  sin2  6 

n  en+l  ’  sin20  \  6n  ’  sin20 


(perpendicular 
polarization)  , 


EnJ  6r4'1  -  sin20  -  £' 


(52) 

rn+l,  n  t  I  t  .  , 

en+l  "  sm 

and 

(53) 

eYndn  =  eandn  ejpndn  ^ 

where 

(54) 

k  tan  ®n 

a  %  — , . =======• 

n  ->  \  J  .  r\ 

n+lv 


en  -  sin2  6 


(!  +  andn)  e 


j^n^n 


(parallel 

polarization), 
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C  J  RICHMOND.  TRANS  AN 0  REFL  COEFF  OF  PLANE  LOWLOSS  MULT  ILAYERS 
C  PEPPENDI CULAR  POLARIZATION  (MAXIMUM  NUMUER  OF  LAYERS  =  8) 

_ DIMENSION,  l>(8>  .L<9)  .(,'(81  .<i(b  )  .SR (9).  TO (8  ). 

2  F  ORM AT  (F 1 0.6.F 1 J .6*+ I0.6.F10.6.F10.6.F I0.O.F1G.6.F10.6) 

3  FORMAT  ( I lb. Fib. 9. FIS. 9. F lb.9.F]b.9) 

A  FORMAT  (Flb.7.Flu.6.Flj.6.Flu.6.FlU.6*Flb.7) 

PF.AQ  3.N.  T.F.OF.FM 

READ  2«D(  I  )  >U(2>  «D(3)«U(  A  ).U(b)  .U(6).l>(7)  .D(O) 

READ  2«E ( 1  )  «t(2) <c(3) >t( A ) .£ (b) «t (6).L(7 )  .E (b> 

PEAU2 • TO ( 1 ).TD(2).TD(J).Tu(A).ru(5).TO(6).r0l7).To(B) 

NN  =  N+l 
E(N+1>  =  I. 

T  =  #0 1 740329*1 

DO  =  6.28318b3*C0S(T  )*(()(  ll  +  D  (  2  )  +0  (  3  !  +D  (  A  )  +D  ( b  )  +0  (6 1  +D  (  7  )  +D 1 8  )  ) 
_S  =  SIN(T)*S|N(T) 

3P( I )  =  SOPT (E ( 1 )-S) 

RO  =  (SR( 1 l-COS(T) )/(oR( 1 )+COS(T) ) 

DO  10  I  =  I ,N 

i i  =  1+1 

SH (111  =  aUHT(t( 11  )-b  > 

_ G(J >  *SH  I  1  ) 

io  ru)  =  (sp( fiT-SR(Tn/(bp( 1 1  )j*sp(iT)  " 

A  =  l.-RR 

_ DO  i5_  l  _= _ 1  »N  . .  .  . 

"~j5  A  =  A#(T.-R(I)) 

A  =  l./A 
M  =  FM 

DO  A 9  J  =  1 . M 
u/  =  1 . 180314E+A/F 
-  ®G  f  G(J  )/W 

. ”'CG"  =  'C0S(GG) 

SG  =  SIN(GG) 

AD  =  3.J  Al_b_926»E  ( 1 )*TD(  1  )«|J(  1  >/(tv*SR<  1  )  ) 

XI  =  CG*(1.-AD> 

Y 1  =  -SG*( 1 .-AD) 

_ X2  =_  -RP+CG*  (  1  .  +AD ) 

...  -yg  -  1rr*sg*(  1  ,+AD) 

X3  =  -RR*CG*( 1 .-AD) 

_ Y3  =  RP»SG* ( 1 .-AO) 

XA  =  CG*(1.+A0) 

YA  =  SG* ( 1 . +AD ) 

00  35  I  =  2.NN, 

I F ( 1 -NN ) 25 • HO • bO 
20  U1  =  1 . 

_ U2  -  -R~(N) 

"u3~  a  -R(N>’~ 

UA  =  1  . 

_ VI  =  _0_. _ 

V2~  =  o". 

V  3  =  0. 

V4  -  0. 

GO "TO  30 
25  II  =  1-1 

_ AD  =  3. 141S926«E( 1)*TD( 1 )*D( )  )/(W*SP( I ) )  _ _ 

GG  =  G ( 1 ) /W 


CG 

= 

COS(GG) 

SG 

= 

SIN(GG) 

U1 

= 

CG*  < 1 .-AD) 

VI 

= 

-SG* ( 1 .-AO) 

U2 

= 

-R< 1 I )*CG*( 1 .+A0) 

V2 

= 

— R (II )*SG*( l.+ADJ 

U3 

= 

-R( I 1 )*CG*(1.-AD) 

V3 

= 

R(  I  I >*SG*( 1 .-AD) 

U4 

= 

CG*( 1 .+AD) 

Fig.  3.  Computer  program  for  transmission  and  reflection 

coefficients  of  low-loww  multilayer  for  perpendicular 
polarization. 
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I 


V4  =  SGMI.+AD) 

30  PI  =  X1*UI-Y1*9I+X2*U3-Y2#V3 
0l_  =  Y1*U|+X|*V1+Y2*U3+_X2»V3 
P2  =  X'|*U2-Y|*V2+X2«U4-Y2»V4 
02  =  YI*U2+X[*V2+X2*V4+Y2*U4 
PJ  =  X3*U1-Y3»VI+X4*U3-Y4*V3 
03  =  Y3*UI+X3*Vl+X4*V3+Y4*U3 
P4  =  X3*U2-Y3*V2+X4*U4-Y4*V4 
04  =  Y3*U2+X3*V2+X4*V4+Y4*U4 
X!  =  PI 
X2  =  P2 
X3  =  P3 
X4  =  P4 
Y1  =  01 
_Y2  a  02 
Y3  =  03 


3b  Y4  =  04 


RCR  =  (-X3*X4-Y3*Y4 >/(X4*X4+Y4*Y4 ) 
PC!  =  l-Y3*X4+X3*Y4 )/(X4*X4+Y4*Y4 ) 
RC2  =  RCR*RCR+RCI*RCI 
RC  =  SORT (RC2 ) 

""  TR" = (X1+X2*RCR-Y2*RC I )*A 
T I  =<Y1+Y2*RCR+X2*RCI )*A 
TC2  =  TR*TR+TI*TI 
TC  *"=  SORT  ( T  C2 ) 

PR  =  ATANIRCI/RCR) 

IF(RCR)40.40.45 
40  PR  =3.1415926  +  PR 

4b  PR  =  b7.295779*PR 

XX  =  ATAN(TI/TR) 

. rF(TR)47.47.4e 

47  XX  =  3.1415926  +  XX 

48  F I  PD  =  -57*295779* (XX+DD/W ) 

‘  "“PUNCH  4.  F  .  RCR  .RC  1 .  RC2.RC  .PR 
PUNCH  4.F.TR.TI .TC2.TC.FIPD 

49  F  =  F  +  DF 

50  STOP . * . 

END 


Fig.  3. 
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C  J  R I  CHI'IQNU*  TRANS  AND  REEL  CutFF  OK  PLANE  LOmLOjb  MULT  ILAYERS 
C  PARALLEL  POLAkUATION  (MAXIMUM  NUMdER  OF  LAYERS  =  b> 

U IMENS'l ON  D ( 5 ) <L ( 6 ) .R< 5 ) .G ( b ) . SR <6>.Tp<5) _ ... 

2  FORMAT  (rTu.6  >F1 0.6.F1  0  .(^fYo.6  .F  I  0.6 .  F  io«6~Fi  0. 6  <  F !  0 .6  f 

3  FORMAT  t | lb. Fib. 9. FIS. V. Fib. 0. Fib. 9) 

_  FORMAT  (F1S.7.F1U.6.F1U.6.FI  J.6  .  F 1  0 .6 .  F  1  S.  7  ) 

READ  3.N. T.F .OF.FM 

READ  2.0(1 >«0<2l«D(3).Dt4).Dl5) 

_  REA0_2.E(n.E(2).E<3>.E(4).E<S)  . .  . .  . . 

READ'  2. TD( 1 ) • TD(2) .TDI3) .TDIA ). TD(5) 

NN  =  N+l 

_  ...  EIN+I  I  =.  1  •  .  _  _ 

T  =  .OI74b329*T 

DD  =  6.2tl3lb33*C0S<T)*(0(  I  ) +0  ( 2  > +0  ( 3  > +0  (  A  > +D  (  b )  ) 

S  »  S|N(T  )*3IM  T  I  .....  _  _ 

SR('i  )  s'sOMTIEU  )-S» 

RR  =  (SR(I)-E(l )*COS(T ) )/(SR( 1 )+E( 1 )*COS(T») 

DO  lu  I  =  l.N 
I  I  =  1  +  1 

SRI  1  I  )  =  SURTIEI I  I >-S> 

G(l)  =  6.2831053*0(1  >*SR( I  ) 

I  U  «(||  =  (Ell  )  *SR  ( I  I )-£ (  I  I  )  *SR  (  I  )  ')'/'(£  (  I  ) *SR  <  1  1  )  +E  (  I  I  )*SR  (  I  )  ) 

A  =  I . -RR 
DO  15  1  =  l.N 
*15  A  =  A*  ( 1  ,-R  (  !  !  ) 

A  =  I./A 
M  =  FM 

DO  49  J  =  I  . M 
41  =  1 . 1BJ314E+4/K 
GG  =  G ( 1  1/W 
CG  =  COSIGG ) 

SG  s  SINIGG) 

_  "  AD  =  3.141 5y26*E  (1  >  * T D  C_1  )#D(  i  )/<  Jl*SR  (1)1 
Xl'  =  CG*"(i.-AD  1  " 

Y 1  =  -SG*(1.-AD> 

Xjl  =  -RR*CG*<  1  .  +AD ) 

.  Y2  =  -HR*S&*  (  i  .  + AD  ) 

X3  =  -RR»CG*( 1 .-AD) 

Y3  -  RR*SG*< 1 .-AD) 

. X4  =  CG*(|.+AD) 

Y4  =  SG* ( I .FAD ) 

_  DO  35  I  =  2.NN  _  _  _ _ 

"!f(T-NN)25 .2  6.50 . . 

.20  Ul  =  1  . 

1)2  =  -R(N) 

U3  =  -R  <  N ) 

U4  -  I  . 

VI  =  0. 

. V2  =  0.' . 

V3  =  0. 

V4  =  0. 

.  GO  TO"  30  . 

2b  II  =  1-1 

_AD  =  3.141 5926*E ( 1 )*TD< I )*D ( I )/ <  to#SR^ I ) ) 

. 7  ~GG  =  G  (T)7i 

""  "  CG  =  COS  (GG ) 

"SG  »  SINIGG) 

“~U1  =*  CG*  (  i . -AD  > 

VI  =  -SG*< I .-AD) 

_U2  =  -R<  Tl  )*CG*(1.+AD)  _ _ 

V2  =  '-R  (  I  f)  *SG*  (  1  .  +AU  I 
U3  =  — R < I  I  )*CG*(1.-AU) 

V3  =  R( II )*SG*( 1 .-AD) 

U4  =  CG* ( V. +AD ) 

Fig.  4.  Computer  program  for  transmission  and  reflection 
coefficients  of  low-loss  multilayer  for  parallel 
polarization. 
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V4  =  SGM1  •  +  AD  ) 

30  PI  =  X1*UI-Y1*VI+X2*U3-Y2*V3 
01  =  YI*U1+XI*VI+Y2*U3+X2*V3 
P 2  =  XI  #U2-Y  |  *V2+X2*U4-Y2*V4 
02  =  Yl*U2+XI*V2+X2*V4+Y2*U4 
P3  =  X3*UI-Y3*VI+X4*U3-Y4*V3 
03  =  Y3*U1 +X3*V I +X4*V3+Y4*U3 
P4  =  X3*U2-Y3*V2+X4*U4-Y4*V4 
04  =  Y3*U2+X3*V2+X4*V4+Y4*U4 
XI  =  PI 
X2  =  P2 
X3  =  P3 
X4  =  P4 
Y I  =  01 
Y2  =  02 
Y3  =  03 
3b  Y 4  =  04 

RCR  =  l-X3*X4-Y3*Y4 >/(X4*X4+Y4*Y4 ) 
RC1  =  1-Y3*X4+X3*Y4 1 / ( X4*X4+Y4*Y4 ) 
RC2  =  RCR*HCR+RCI*RCI 
«C  =  SORT ( RC2 ) 

TR  =<X1+X2*RCR-Y2*RCI )*A 
T1  =  (Y1+Y2*RCR+X2*RCI )  *A 
TC2  =  TR*TR+TI*TI 
TC  =  SORT  t  TC2 ) 

PR  =  ATANiRCI/RCR) 

IF<RCR>40,40,45 
40  PR  =  3.1415926  +  PR 

45  PR  =  57.295779*PR 

XX  =  ATANITI/TR) 

I  FUR  147*  47. 48 
47  XX  =  3.1415926  +  XX 
40.  FI  PD  =  -57.295779* I XX+DD/W 1 

PUNCH  4.F.HCR.RCI .RC2.RC.PR  " 

PUNCH  4.F.TR.n,TC2.TC.FlPD 

49  F_=  F_+  OF _ 

50  STOP  “  '  . * 

END 


Fig.  4. 
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and 


(55)  pn  *  k  Jt#n  -  sin20  . 

In  Eqs.  (51)  through  (55),  tn  and  represent  the  real  part  of  the 
relative  complex  permittivity  of  layers  n  and  n+1.  The  errors  in  the 
transmission-coefficient  data  arising  from  these  approximations  are  of 
second  order,  and  the  resulting  data  are  accurate  for  low-loss  multilayers 
such  as  radome  sandwiches. 

The  symbols  used  in  the  computer  programs  in  Figs.  3  and  4  are 
defined  in  Table  I. 


TABLE  I 

Symbols  Used  in  the  Computer  Programs  in  Figs.  3  and  4 

N  =  number  of  layers 

T  =6,  angle  of  incidence,  degrees 

F  =  lowest  frequency,  megacycles 

DF  =  increment  in  frequency,  megacycles 

FM  =  number  of  frequencies  to  be  calculated 

D( I)  =  thickness  of  layer  I,  inches 

E(I)  =  real  part  of  relative  permittivity  of  layer  I 

TD(I)=  loss  tangent  of  layer  I 

S  =  sin26 

SR(I)  =  Je'T  -  sin2  6 

G(I)  =  Ztt  djJ e j  -  sin20 

R(I)  =  interface  reflection  coefficient,  rj^j 
W  =  wavelength  in  free  space,  inches 

The  quantities  calculated  and  punched  as  output  data  are  given 
in  Table  II. 

These  computer  programs  were  used  with  an  IBM  1620  digital 
computer  to  study  the  broadband  characteristics  of  a  sandwich  with  five 
layers  having  dielectric  constants  of  5,  4,  3,  2,  and  1.5  and  thicknesses 
of  0.15,  0.15,  0.15,  0.26,  and  0 . 26  inches .  The  results  for  frequencies 
from  zero  to  40  gigacycles  are  shown  in  Figs.  5,  6,  and  7  for  the  lossless 
and  the  low-loss  cases.  Each  of  the  six  curves  involved  calculation  of  the 
transmission  coefficient  at  80  frequencies.  The  calculation  time  was 
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Fig.  5.  Power  transmission  coefficient  versus  frequency 
for  five-layer  sandwich. 


Fig.  6.  Power  transmission  coefficient  versus  frequency 
for  five-layer  sandwich. 
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Fig.  7.  Power  transmission  coefficient  versus  frequency 
for  five -layer  sandwich. 


15  seconds  for  each  frequency,  which  amounts  to  20  minutes  for  each 
curve.  It  is  apparent  from  Fig.  6  that  this  particular  sandwich  has 
excellent  broadband  properties  for  parallel  polarization  at  oblique 
incidence . 


T7TBLE  II 

Symbols  for  Output  Data  of  Computer  Programs 
in  Figs •  3  and  4 

R  =  RCR  +  j  RCI  =  reflection  coefficient 
RC2  =  f  R  | 2  =  power  reduction  coefficient 
RC  =  [  R f  =  voltage  reflection  coefficient 
PR  =  reflection  pha:  angle,  degrees 
Tn  =  TR  t  j  TI  -  normal  transmission  coefficient 
TC2  =  f  Tn  | 2  =  power  transmission  coefficient 
TC  =  |  Tn  |  =  voltage  transmission  coefficient 
FIPD  =  insertion  phase  delay,  degrees 
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When  used  with  the  IBM  1620  equipment  at  The  Ohio  State  University, 
these  programs  handle  a  maximum  of  eight  layers  for  perpendicular  polari¬ 
zation  and  five  layers  for  parallel  polarization.  With  slight  modifications 
in  these  programs,  the  maximum  number  of  layers  can  be  increased  by 
using  equipment  with  greater  storage  capacity.  If  the  reflection  phase,  the 
transmission  phase,  and  the  insertion  phase  delay  need  not  be  calculated, 
the  maximum  number  of  layers  can  be  increased,  even  while  using  the 
IBM  1620  equipment,  by  deleting  the  statements  which  call  for  the  arc¬ 
tangent  subroutine. 


III.  SURFACE  WAVES  ON  LOSSLESS 
PLANE  MULTILAYERS 

Consider  a  surface  wave  propagating  on  a  plane,  lossless  multi¬ 
layer.  Let  h  represent  the  phase  constant  of  the  surface  wave,  and  let 
a  represent  the  normal  attenuation  constant  in  the  free-space  regions 
on  both  sides  of  the  multilayer.  Both  h  and  Q  are  real  quantities.  For 
a  TE  surface  wave,  the  field  in  the  free-space  regions  on  both  sides  of 
the  multilayer  has  the  form 


(56) 


Ex  =  e-Jhye-a,z| 


The  multilayer  and  the  coordinate  system  are  shown  in  Figs.  1  and  2. 
The  electric  field  intensity  Ex  in  layer  n  is  given  by 


,,7.  „  ,  -^n*  V.  -Jhy 

(57)  En  =  (a  e  +  b  e  )  e 


where 


(58) 


Yn  =  kjh2/k2  -  pnen  . 


The  quantity  under  the  square-root  sign  in  Eq.  (58)  may  be  positive  in 
some  layers  and  negative  in  others.  Thus,  the  field  in  some  layers  is 
"evanescent"  while  in  other  layers  it  is  the  sum  of  two  criss-crossing 
plane-wave  fields. 
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It  can  be  shown  that 


(59)  h/k  =  \/\g  =  c/v 

where  X  is  the  wavelength  in  free  space,  c  is  the  velocity  of  light  in 
free  space,  \g  is  the  "guide  wavelength,  "  and  v  is  the  phase  velocity  of 
the  surface  wave. 

The  equations  for  a  surface  wave  on  a  plane  multilayer  differ  from 
those  in  Section  II  only  in  that  the  quantity  k  sin  0  is  replaced  by  h,  and 
in  the  surface-wave  case  the  amplitude  of  the  incident  plane  wave  is  zero. 
Thus,  the  matrix  Eq.  (25)  applies  also  for  surface  waves,  for  which  case 

(60)  D0  =  0 
and 


(61)  An+1  =  0  . 

Now  Eq.  (25)  becomes 


from  which 

(63)  dBN+i  =  0. 

But  B^-fi  represents  the  amplitude  of  the  surface-wave  field  at  one  of 
the  outer  surfaces  of  the  multilayer,  and  it  may  be  assumed  to  be  non¬ 
zero.  Hence, 

(64)  d  =  0  . 

Equation  (64)  represents  the  "transcendental  equation  for  surface 
waves,  "  since  d  represents  the  lower  right  element  in  the  matrix 
(1/t)  Mi  •  M2  •  •  •  Mjq  ’  S  (see  Eq.  (25))and  the  equation  d  =  0  cannot  be 
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rearranged  to  obta'n  an  explicit  expression  for  h.  It  is,  however,  a 
useful  implicit  equation  for  h. 

Since  the  zeros  of  the  factor  1/t  are  not  signuicant  here,  the 
transcendental  equation  for  surface  waves  can  be  written  as  follows: 


(65) 


Mj 


•  M2  •  M3  •  •  •  •  S  = 


b' 

d 


in  which  d  is  required  to  vanish;  and 
/  "8n  ..  8n' 


(66) 


Mn  = 


-gn  8n 

l*  a  ^ 11  O 


■  -rne 


(67) 


S  =1 


.  -r 


N+l 


■rN+r 
i  j 


(68) 

and 


gn  ~  ^n  ^n  > 


(69) 


H’n-l  Yn  "  M-n  m-1 


rn  =  rn,  n-1-  v  ,  u  v  ,  (for  TE  waves) . 
r  n- 1  *n  f  'n-l 


It  can  be  shown  that  the  range  of  h/k  is  from  unity  to  J  (JLmem  where 
pm  and  em  are  the  relative  permeability  and  permittivity  of  the  layer 
having  the  greatest  value  of  the  product  pe.  If  h/k  =  1,  the  surface  wave 
is  said  to  be  at  cutoff.  The  value  h/k  pmem  is  the  high-frequency  limit. 
To  determine  the  phase  constants  of  all  the  surface  waves  that  may  propa¬ 
gate  on  a  given  multilayer,  one  may  calculate  the  lower-right  matrix 
element  d/  in  Eq.  (65)  for  many  values  of  h/k  in  the  range  from  unity  to 
J  pmern'  The  zeros  of  d*  can  then  be  found  by  making  a  graph  of  |  d;  | 
versus  h/k,  or  by  finding  the  roots  of  a  polynomial  which  fits  the  calculated 
values  of  d7. 
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The  matrix  multiplications  indicated  in  Eq.  (65)  arc  programmed 
in  Fortran  as  shown  in  Fig.  8.  The  symbols  in  this  program  are  defined 
in  Table  III.  The  program  applies  to  nonmagnetic  multilayers  (|i  =  p0  an<^ 
H-r  =  l>- 


TABLE  III 

Symbols  Used  in  the  Computer  Program  in  Fig.  8 

N  =  number  of  layers 

M  =  number  of  values  of  h/k  to  be  used 

T  =  initial  value  of  h/k,  on  input  data  card 

T  =  current  value  of  h/k,  in  output  data 

DT  =  increment  in  h/k 

F  =  frequency,  me 

D(I)  =  thickness  of  layer  I,  inches 

E(I)  =  relative  dielectric  constant  of  layer  I 

W  =  wavelength  in  free  space,  inches 

T2  =  (h/k)2 

d'  =  X4  +  j  Y4 

Z4  =  |d'|  . 

Figure  9  is  a  graph  of  the  data  obtained  with  this  program  on  an 
IBM  1620  digital  computer  for  a  sandwich  having  three  layers  with 
dielectric  constants  of  4,  1,  and  4,  and  thicknesses  of  0.25,  0.20, 
and  0.25  inches.  From  Fig.  9,  this  sandwich  supports  TE  modes 
having  h/k  =  1.62  (TEi  mode)  and  1.70  (TE0  mode)  at  11,803  me.  This 
agrees  accurately  with  independent  calculations. 

This  computer  program  was  also  used  to  obtain  surface-wave 
data  for  the  five-layer  sandwich  described  in  Figs.  5,  6,  and  7.  At 
10  gigacycles  the  only  TE  modes  it  supports  are  the  TE0  and  TEi  modes, 
for  which  h/k  is  1.85  and  1.30,  respectively. 

Since 


(70)  d  -  t0,i  til2t2, 3  •••  tj\ji  n+i  d 

and 

(71)  In,  n-f-1  _  1  +  rn,  n+1  *  rn+l,  n 

[d  |  has  some  zeros  which  |d|  does  not  have.  These  extraneous  zeros 
occur  at  values  of  h/k  for  which  one  of  the  interface  reflection  coeffi¬ 
cients  rn  takes  on  the  value  +1.  From  Eqs.  (58)  and  (69),  this  occurs 


f -i 
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C  J  R I CHMOND.  SURFACE  A'AVES  ON  PLANE  LOSSLESS  MULT  ILAYER.  TE 

C MAX  I  MUM  NUMBER*  OF  LAYERS .  »  .  U _ _ _ _ 

dimension  DCBiTi'c'ai 

2  FORMAT  (FlJ.6*Flu.6*FlL.b*F10.6.K|o.6iF10.6.FI0.6*F|0.6) 

3  FORMAT  (I  15.  I  l5.FJ6.9iFI5.9.F|5.y ) _ _ 

4  FORMAT  lFI0.V.Fli.9.H‘J.y<F|^.y.FIS.y) 

WtAD  3.N.M. T.OT.F 

READ  3. 01ilt0(2!iD(3)<D(AI. 0(d). DI6). 0(7)10(0) 

READ  2.E(l).E(£).E(J).EI«).c(b).t(6).E(7).E(0) 

N I  =  N+l 

PI  =  3.14  16926  .  . . . 

T  =  T-L)T 

V  =  I . 16l3)4£+4/F 
DO  70  L  =  I.M 

T  =  T  +  or 
T  2  =  T*T 

SS  =  SORT ( T2- 1 .  ) 

A  =  AiJS(T2-L  (1)1 
S  =  SORT ( A ) 

G  =  2.*P|*U( 1 )*S/W 
IF ( T2-E1 1  l )  10. 10.  1  I 

10  GR  =  COS (G ) 

G 1  =  S I N  (G ) 

FR  =  GR 

FI  =  -G I 
GO  TO  15 

11  GR  =  EXP(G) 

G I  =  0. 

FR  =  1 ./GR 

FI  =  0. 

15  IF(T2-E<1 1)16.16.17 

16  PH  =  Pl-2.*ATAN<5/3s> 

RR  =  COS (PH) 

R 1  =  SIN(PH) 

.  GO  TO  20  .  .  ....  _ 

17  RR  =  (S-SS)/(S+SS> 

R I  =  0. 

20  XI  =  FR . . _  _ _ _ 

Y 1  =  FI 

X2  =  -RR*GR+RI+GI 

.  Y2  =  -R I  «GR-RH*G  I  . . .  . . . . 

X3  =  -RR*FR+RI»FI 

Y  3  =  -R I *FR-RR*F 1 
X4.  =  GR_ 

"Y4~""G] 

DO  65  1  =  2.N1 

_ I F ( 1 —Ml  135.25.86 _  _ 

25  A  =  Af)S(T2-E(N)  ) 

S  =  SORT (A) 

IF  (T2-E1N)  )26.26.J!7 _ _  _ _ 

26  PH  =  2.*ATAN(S/SS) 

RR  =  COS (PH) 

R I  =  -SIN(PH) 

GO  TO  '30  ~  “ 

27  RR  =  (SS-S)/(SS+S) 

Rj  J=_U_.  _ 

30  ui  =  j."_’ 
vi'  =  0. 

U2  =  -RR 
V2  =  -R 1 
U3  =  -RR 

V3  =  -R I _ 

U4  =  1  . 


Fig.  8.  Computer  program  lor  phase  constants  of  TE 
surface  waves  on  plane  multilayers. 
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3b 


VA  -  u. 

GO  TO  60 

A  =  AUSI T2-E (1)1 
)S  =  SORT (A) 

,  G  =  2.  »PI  »0(  I  )»5/*' 

'  II  =  l-l 
A]  =  A3SIT2-EI  ID) 

SI  =  SORTIAI  ) 
i  -  IT2-EI I  1 )»<T2-E(  III! 
IFI2lAO.36.36 
36  RR  =  (S-Sl 1/IS+S1 1 
R I  =  0. 

GO  TO  bO 


AO  I F I T2-E (I11A1.AI.A5 


A  1 

Pm 

=  PI-2.*ATANIS/Sl 1 

RR 

=  COS (PH) 

R  1 

=  SINIPH) 

GO 

TO  50 

4b 

PH 

=  2.*ATAN(SlXS) 

RR 

=  COSIPH) 

R  1 

=  -SINIPH) 

50 

IFIT2-EI 1 1 151.51 .52 

51 

GR 

=  COS(G) 

G 1 

=  SINIG) 

FR 

=  GR 

F  1 

=  -G  1 

GO 

TO  55 

52 

GR 

=  EXPIG) 

G I 

=  0. 

FR 

=  l./GR 

FI 

=  J. 

55 

01 

=  FR 

VI 

=  FI 

U2 

=  -RW#oR+R!*Gl 

V2 

=  -RI*GR-RR*GI 

U3 

=  -RR*FR+R1#FI 

V3 

=  -RI*rR-RR*FI 

U4 

=  GR 

V4 

=  G I 

60 

P! 

=  X|*U1-Y1*VI+X2*UJ-Y2*V3 

01 

=  Y1*01+X|*V1+Y2*U3+X2*V3 

P2 

=  X1*U2-Y1#V2+X2*U4-Y2*V4 

02 

=  Y1*U2+X1*V2+Y2*UA+X2*VA 

P3 

=  X3*U1-Y3*V1+XA*U3-YA*V3 

03 

=  Y3*U1+X3*V1+YA*U3+XA*V3 

P4 

=  X3*U2-Y3*V2+XA*UA-YA*V4 

Q4 

=  Y3*U2+X3*V2+YA*U4+XA*V4 

XI 

=  PI 

X2 

=  P2 

X3 

=  P3 

X4 

=  P4 

VI 

=  01 

V? 

=  02 

Y3 

=  03 

65 

Y4 

=  QA 

ZA 

=  SORT <XA»XA  +  YA*YA  1 

70 

PUNCH  4»T»X4»V4»2:4 

30  STOP 
END 


Fig.  8. 


1751-2 


23 


Fig.  9.  Matrix  el 
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d'  versus  h/k  for  three-layer  sandwich. 


when  Yn  =  0  or 

(72)  h/k  =  J|in6n  • 

These  zeros  of  Id^  are  not  significant  with  respect  to  surface  waves. 
They  are  easily  identified,  thus  permitting  the  significant  zeros  to  be 
recognized.  For  example,  in  Fig.  9  the  zeros  at  h/k  =  1  and  2  arc  to 
be  ignored. 

TM  surface  waves  on  a  plane  multilayer  can  be  analyzed  by  means 
of  a  similar  program  using  the  expression  for  the  interface  reflection 
coefficient  for  parallel  polarization  instead  of  Eq.  (69). 

IV.  TRANSMISSION  COEFFICIENT  OF 
AN  INHOMOGENEOUS  LAYER 

Consider  a  plane  wave  incident  on  a  plane  inhomogeneous  layer 
immersed  in  free  space,  as  in  Fig.  10.  If  the  permittivity  is  a  function 
of  the  normal  coordinate  only 

(73)  e=e(z) 
and 

(74)  H=Po- 
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Using  step-by-step  numerical  integration,  equations  are  developed 
in  Report  1180-44  for  the  electric  field  intensity  distribution  Ex(z)  and 
the  transmission  coefficient  of  such  a  layer.  These  equations  are  listed 
below  for  perpendicular  polarization. 

If  the  electric  field  intensity  is  set  equal  to  unity  (E0  =  1)  at  the 
coordinate  origin,  the  field  at  the  point  (x,y,z)  =  (0,0,  h)  is  given  by 


(75)  Ei  =  1  -  k2h2(e,  -  sin2 0) / 2 

T 

T  i-2  h  2  k3h3 

+  j  j^kh  cos  0  +  ——  e+  tan  6+  -  ~j~  (e+  -  sin  0)  cos  0 

where  e^.  and  tan  6+  represent  the  relative  permittivity  and  the  loss  tangent 
just  inside  the  inhomogeneous  layer  at  z  =  0+. 


Let  En_i,  En,  and  En^q  represent  the  electric  field  intensities  Ex 
at  the  points  z  =  (n-l)h,  nh,  and  (n+l)h.  Further,  let  R  and  S  represent 
the  real  and  imaginary  parts  of  the  electric  field  intensity;  for  example, 

(76)  En  =  Rn+jSn. 

The  recursion  formulas  for  the  real  and  imaginary  parts  of  the  electric 
field  intensity  in  the  inhomogeneous  layer  are 


(7?)  Rn+1  =  t2  -  k2h2(en  -  sin20)]Rn  -  k2h2enSn  tan  6n-Rn.1 

and 

(78)  Sn+1  =  [Z  -  k2h2(en  -  sin20)JSn  +  k2h2enRntan  5n  -  Sn_1 

where  en  and  tan  6n  represent  the  relative  dielectric  constant  and  the 
loss  tangent  at  the  point  (x,  y,  z)  =  (0,  0,  nh) . 

After  Eq.  (7  5)  is  used  to  calculate  the  field  Ei  at  z  =  h,  Eqs.  (77) 
and  (78)  may  be  employed  to  calculate  the  field  (R2  and  S2  )  at  z  =  2h. 
Following  this,  Eqs.  (77)  and  (78)  are  used  again  to  calculate  the  field 
at  the  next  point  z  =  3h,  and  so  on  until  the  inhomogeneous  layer  has 
been  traversed. 
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Finally,  the  transmission  coefficient  is  determined  by  the  following 
expression: 


(79) 


4  kh  cos  0 

kh  cos  0  -  j(Ejsj+j  -  E^j_j) 


where  Ej^j,  and  Ejyj+|  are  the  electric  field  intensities  at  the  points 

z  =  d-h,  z  =  d,  and  z  -  d+h,  with  d  representing  the  thickness  of  the 
inhomogeneous  layer. 


The  step  size  h  must  be  small  to  obtain  accurate  results  with 
these  equations. 


Figure  11  shows  a  computer  program  based  on  the  above  equations, 
table  IV  defines  the  symbols  used  in  this  program. 

TABLE  IV 

Symbols  Used  in  Computer  Program  in  Fig.  11 

N  =  number  of  field  points  to  be  calculated 
H  =  kh 

TH  =  increment  in  0,  degrees 
THETA  =  angle  of  incidence  6,  radians 
K  =  angle  of  incidence  0,  degrees 
S  =  sin20 
C  =  cos  0 
H2  =  (kh)2 

E(I)  =  relative  dielectric  constant  at  point  I 
R1,R2,R3  =  real  parts  of  electric  field  intensity 
at  points  1,  2,  and  3 

S1,S2,S3  =  imaginary  parts  of  electric  field  intensity 
at  points  1,  2,  and  3 

EM  =  |E| 

TR  +  j  TI  =  complex  transmission  coefficient 
TS  =  power  transmission  coefficient 
T  =  voltage  transmission  coefficient  magnitude 
PH  =  phase  of  normal  transmission  coefficient 


Some  of  the  results  obtained  with  this  digital  computer  program 
are  listed  in  Tables  V  through  VIII.  These  results  are  found  to  agree 
very  well  with  the  exact  solutions  for  homogeneous  and  exponentially  - 
inhomogeneous  layers. 
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C  J  RICHMOND.  LoSSLESo  1  NriyH  JOt  NiOUS  PLANt.  LAVcR.  PerPEND I  GULAK 

~C  I'LXCWC  FI  C'lO  DI  STRIUDT  ION  A  No  TRANS-'!  I  SS I  ON  "COEFFICIENT 
DIMENSION  i  ( 1 6  ) 

2  FORMAT  <  l3.K|5.U.F15.tl.KI3.b.FI5.U.F15.H> 

3  FORMAT  (F10.6.F|9.6»Flij.o.Fl'J.6.FlJ.6.l*  ID.6.FIG.O.F  10.6) 

RCAO  2.N.M.TH 

H2  -  H»H 

READ  3  1  )'."E<2>  .A  (3)  W..I4  )  .r.Cbl  tE(6  )  .r.(7)  .t.  16) 

READ  3.L(9>.t<lu).ECIl  ) • l ( 1 2  ) . h 11 3 ) • t: II  A ) . L  il 6 ) t  l (  I  6 ) 

M  =  9C./TH 

TH"  =  .Oi  7o'?)J29+TH 
THETA  =  0. 

DO  26  J  =  Ii M 

. S"'=  S  i  N(  THETA  j*S  INI  TriETA) 

C  =  COS ( THETA  1 

R1  =  1  . 

51  -  9 . 

I  =  1 

PUNCH  2  < I *R 1 «SI .R) 

A  *'ETn  -"3 

R2  =  1.  -  N2»A/2.  -  H2*(c(2)-l.(  I  )  I/O. 

52  =  Hi«C-H2*H*A*C/6. 

I  =~2  . 

EM  =  SORT (R2*R2+52*S2 ) 

PUNCH  2  •  I  . R2 . 32 . EM 
-  DO'TO  I  "5  2.N- 

II  =  1+1 

A  =  Ell  )  -3 

""3  •S"r.-K2»A 
R3  =  3*W2-R| 

53  =  6*32-31 

-  EM"  TF  •S0RTrR3»R!j>S3*S3 )  . . — . 

PUNCH  2. II .R3.S3.EM 
I F { N- I >30.15.5 

-  SI  A  S3 

R2  =  R3 

TO"  S2"=“53" . . .  . . . . 

15  DR‘  -  2.*H*C*R2+S3-S1 
D I  =  2.*H*C*S2-R3+R1 

. .  TR  =  4S*H*C*T)R/IDR*DR+DI*DI  ) 

T I  =  -4.*H*C*DI/(DR*f;R  +  JI*Dl  ) 

TS  =  TR*TR+TI*TI 

. . . T  "  =""50RT "( T5T  " . ~  "  . 

PH  =  ATANITI/TR) 

[FITR120.25.25 

- 2TT"PFrs— PH+T.TA  1 3TOET .  . . .  . 

.  25  PH  =‘ PH*B7. 29570 

"  K  =  57 .2930* THETA 

. -P0n£FT2'VK  t  TR'.T  f .  T'S.  T  .PH"  ' 

punch  2 

26  THETA  =  THETA  +  TH 
FND 


Fig.  11.  Computer  program  for  field  distribution  and  transmission 
coefficient  of  inhomogeneous  plan^ layer  for  perpen¬ 
dicular  polarization. 
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TABLE  V 


Electric 

Field  Intensity  Distribution 

in 

Homogeneous 

Layer 

|E 

1 

kz 

ix 

0  =  0 

0  =80° 

0 

4 

1.0000 

1.0000 

.  1 

4 

.9850 

.9850 

.2 

4 

.9411 

.9405 

-3 

4 

.8717 

.8678 

.4 

4 

.7827 

.7693 

.5 

4 

.  6835 

.6482 

.  6 

4 

.5888 

.5086 

.7 

4 

.5199 

.3561 

.8 

_ 

.4999 

.2011 

TABLE  VI 

Electric  Field  Intensity  Distribution 
in  Inhomogeneous  Layer 


e_  =  4  exp(-. 3077  kz) 

_ jEj _ 


kz 

£r 

O 

II 

CD 

0  =  80° 

0 

4.0000 

1.0000 

1.0000 

.  1 

3.8787 

.9852 

.9852 

.2 

3.7612 

.9427 

■9?20 

.3 

3.6473 

.8768 

.8730 

.4 

3.  5367 

.7938 

.7810 

.  5 

3.4296 

.7030 

.  6696 

.  6 

3.  3257 

.6171 

.5428 

.7 

3.  2249 

.5533 

.4052 

.8 

— 

.5297 

.2636 
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TABLE  VII 


y 

Voltage  Transmission  Coefficient  of  Homogeneous  Layer 
for  Perpendicular  Polarization 
er  =  4,  d  =  .111408  \Q 

I  T  | 


_0_ 

Computed 

Exact 

0° 

.8076 

.8042 

10 

.8023 

.7990 

20 

.7859 

.7827 

30 

.7560 

.7530 

40 

.7087 

.7060 

50 

.6379 

.6356 

60 

.  5354 

.5335 

70 

.3934 

.3921 

80 

.2106 

.2099 

TABLE  VIII 

Voltage  Transmission  Coefficient  of  Inhomogeneous  Layer 
for  Perpendicular  Polarization 
er  =  4  exp(-,3077  kz) 
d  =  .111408  \0 


9  T  Computed 


0° 

.8357 

10 

.8310 

20 

.8161 

30 

.7889 

40 

.7450 

50 

.6774 

60 

.5761 

70 

.4298 

80 

.2331 

While  the  calculated  results  listed  in  these  tables  are  sufficiently 
accurate  for  most  radome  design  purposes,  more  accurate  data  can  be 
obtained  quite  readily  by  using  the  program  shown  in  Fig.  11  with  a 
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smaller  value  of  the  step  size  h  or  kh.  With  an  IBM  1620  digital  com¬ 
puter,  a  total  running  Mine  of  six  minutes  was  required  to  obtain  the 
data  listed  in  Tables  V  t  rough  VIII. 


V.  CONCLUSIONS 

Digital  computer  programs  are  given  for  the  transmission  and 
reflection  coefficients  of  low-loss,  multilayer,  radome  sandwiches. 
Using  matrix  multiplication,  these  programs  yield  design  data  for 
parallel  and  perpendlc.  ,ar  polarization  in  the  most  efficient  manner. 
Typical  data  are  included  for  a  five -layer  sandwich  for  normal  and 
oblique  incidence  at  frequencies  from  zero  to  40  gigacycles. 

Matrix  multiplication  is  also  employed  to  calculate  the  phase 
velocity  of  the  surface  waves  which  may  propagate  on  a  plane  multilayer 
sandwich.  Typical  results  are  given  for  three-layer  and  five-layer 
sandwiches . 

The  transmission  coefficient  of  an  inhomogeneous  plane  layer  is 
determined  with  step-by-step  numerical  integration  of  the  wave  equation. 
A  computer  program  is  included  for  this  problem,  along  with  typical 
computed  data  for  an  exponentially-inhomogeneous  layer. 
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